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ABSTRACT 

Aims. We performed a detailed study of maser and radio continuum emission toward the high-mass star-forming region G23. 01-0.41. 
This study aims at improving our knowledge of the high-mass star-forming process by comparing the gas kinematics near a newly 
bom young stellar object (YSO), analyzed through high spatial resolution maser data, with the large-scale environment of its native 
hot molecular core (HMC), identified in previous interferometric observations of thermal continuum and molecular lines. 
Methods. Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful 
maser species in G23.01-0.41: HjO at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). 
In addition, we performed high-resolution (> O'.'I), high-sensitivity (< 0.1 mjy) VLA observations of the radio continuum 
emission from the HMC at 1.3 and 3.6 cm. 

Results. We have detected H2O, CH3OH, and OH maser emission clustered within 2000 AU from the center of a flattened HMC, 
oriented SE-NW, from which emerges a massive '^CO outflow, elongated NE-SW, extended up to the pc-scale. Although the three 
maser species show a clearly different spatial and velocity distribution and sample distinct environments around the massive YSO, 
the spatial symmetry and velocity field of each maser specie can be explained in terms of expansion from a common center, which 
possibly denotes the position of the YSO driving the maser motion. Water masers trace both a fast shock (up to 50 km s"' ) closer 
to the YSO, powered by a wide-angle wind, and a slower (20 km s"' ) bipolar jet, at the base of the large-scale outflow. Since the 
compact free-free emission is found offset from the putative location of the YSO along a direction consistent with that of the maser 
jet axis, we interpret the radio continuum in terms of a thermal jet. The velocity field of methanol masers can be explained in terms 
of a composition of slow (4 km s"' in amplitude) motions of radial expansion and rotation about an axis approximately parallel to 
the maser jet. Finally, the distribution of line of sight velocities of the hydroxyl masers suggests that they can trace gas less dense 
(uh, < 10' cm"') and more distant from the YSO than that traced by the water and methanol masers, which is expanding toward the 
observer. A few pairs of OH masers, with different circular polarization, are well aligned in position on the sky and we interpret them 
as Zeeman pairs. From Zeeman splitting, the derived typical values of the magnetic field are of a few mG. 

Key words. Masers - Techniques: high angular resolution - ISM: kinematics and dynamics - Stars: formation - Stars: individual: 
G23.01-0.41 
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1 . Introduction statistically found far away from the observer, at distances of a 

• few kpc, clustered in tight associations. 
T— I . Hot, dense, molecular cores (HMCs) of dust and gas located 

L; ■ in giant molecular clouds (GMCs) are the birth sites of high- Since a few years, we have started an observational cam- 

."^ ' mass young stellar objects (YSOs). Ultraviolet light from the P^ign to study the high-mass star-forming process by compai-- 

X newly formed stars ionizes the surrounding gas, creating a com- ing interferometric thermal date, tracing the l^ge-scale environ- 

U ■ pact H n region (e.g.,jHoare et al. 2007), excites several molecu- ment (e.g., | Codella et al. || 1997t LFuruya et al."2008'), with Very 

."^ ; lar transitions (e.g., ISridharan et al.n20021 . and heats the dust. Long Baseline Interferometry (VLBI) measurements of maser 

whic h reradiates the ener gy in the far infrai-ed (FIR) band transitions, tracing the inner kinematics of the (proto)stellar co- 

(e.g., ' Molinari et al.l Ho^l). At the radio wavelengths, maser coon. Deteils about our VLBI observational program to mea- 

emissions of several molecular transitions are observed during sure molecular masers in HMSFRs are extensively presented in 

the earliest evolutiona ry phases of high-mass (proto)stai-s (e.g., ISanna et al. ('2010, hereafter Paper I). The present paper focuses 

ISzymczaketall l2Uol . even before the appeai-ance of an ul- on our observations and analysis of the HMSFR G23.01-0.41. 

tra compact H 11 region (UCH 11). From an observational point In Sect. 2, we provide an up-to-date review of previous ob- 

of view, however, the study of high-mass star-forming regions servations toward the HMSFR G23.01-0.41. Section 3 describes 

(HMSFRs) is challenging because of three main limitations: the our VLBI observations of the 22.2 GHz water (H2O), 6.7 GHz 

dusty environment is obscured at optical/NIR frequencies; time methanol (CH3OH) and 1.665 GHz hydroxyl (OH) maser tran- 

scales of massive star formation ai-e short, and hence, the chance sitions, together with the new Very Large Array (VLA) obser- 

of observing massive YSOs is small; massive YSOs ai-e rare and vations of the radio continuum emission at 1.3 and 3.6 cm. In 

Sect. 4, we illustrate the spatial morphology, kinematics, and 

Send offprint requests to: A. Sanna, e-mail: asanna@ca . astro . it time-variability of individual maser species, and present results 



2 



A. Sanna et al.: VLBI study of maser kinematics in high-mass SFRs. II. G23. 01-0.41 



from our VLA observations, constraining the properties of the 
radio continuum observed associated with the masers. Section 5 
discusses the spatial association of the maser species and their 
overall kinematics, and draws a comprehensive picture of the 
phenomena observed in the HMSFR G23.01-0.41 on angular 
scales from a few mas to tens of arcsec. Main conclusions are 
summarized in Sect. 6. 



2. The HMSFR G23.01-0.41 

The HMSFR G23.01-0.41 is loc ated at a distance of 
4. 59^0 33 kpc dBrunthaler et al.ll2009l) . has a large-scale clump 
mass of about 4x10"^ Mq 

a5 



('Furuv a et al.ll2008h . and a bolometric 
luminosity o f about 10^^ (upper limit) inferred from IR data 
dArava et all 12008 fl. From observations of NH3 and CH3CN 
transitions, the systemic velocity of the region (Vsys) with re- 
spect to the local stan dard of rest (LSR) was determ ined to be 
77.4kms-i (.Codella etan[T997l: iFuruva et alJl2008h . 

On an arcminute scale, Spitzer IRAC GLIMPSE observa- 
tions show strong 4.5 fim excess indicative of shocked gas 
(lArava et al.l l2008l and references therein), at the position of 
a molecular clump extended over an area of about 30" x 
30", imaged with the Nobeyama Millimeter A rray (NMA) 
in the '^CO (1 - 0) and C"^0 (1 - 0) lines (iFumva et alJ 
l2008h . Associated with the extended mid-IR emission, a mas- 
sive '^CO (1-0) bipolar outflow was observed, with a promi- 
nent red-wing emission, elongated in the NE-SW direction. 
The same outflow was also det ected in the ^^CO (1-0) and 
HNCO (5o5 - 4o4) lines (Furuva et alj|2008l) . Single-dish ob- 
servations toward G23. 01-0.41 detected several, bo th thermal 
and maser, molecular transitions. 'Anglada et al.l (1199 6) reported 
similar spatial and velocity distribution of the CS (1 - 0) 
and NH3 (1,1) lines tracing the turbulen t, high-density (10** - 
10^ cm"^) gas in the region (see also iLarionov et al] 1 19991) . 
Thermal, broad line (FWHM> 20 km s~' ) SiO emission in 
the ground vibrational stat e (J = 2 - 1 and 3 - 2) was de- 
tected by Harju et al. ( 1998). ICaswell et al.l(l2000l) d etected ther- 
mal C H3OH emission from the 156.6 GHz (see also lSlvsh et alj 
Il999l) and 107.0 GHz transitions, this latter superimposed on a 
maser compone nt. Methanol maser emission was al so reported 
at 6. 7 GHz (e.g.,|M enten 1991; Goedhart et al."2004'), 12.2 GHz 



(e.g.JCasweU et aL1993; Blaszkiewicz & Kus 2004^44.1 GHz 



(ISlvshet al.l ll994r and 95.2 GHz (TVal'tts et al.li^ 



First at- 
tempts to measure Zeeman splitting of the 6.7 GHz maser com- 
ponents suggest a magnetic field strength of a few tens of mG 
(IVlemmingsll2008h . The G23. 01-0.41 re gion hosts both strong 
22.2 GHz H2O (e.g. JCasweU et alJll983tls"zvmczak et al.ll2005n 
and 1.665, 1.66 7 (main-lines), and 1.612 GHz (satellite-line) 
OH masers (e.g. .'Caswell & Ha vnes|[T983l : ISzvmczak & GerardI 
[2004). Throu gh Z eeman splitti ng of the 1.667 GHz maser com- 
ponents, !Szjanczak&_Gerard| (120041) estimated a magnetic field 
strength of a few tenths of mG. Interestingly, in the spec- 
tra of the 1.665, 1.667, and 1.720 GHz OH transitions, ab- 
sorption features were also noted ("Caswell &Havned 1 1983t 
[Szymczak & Gerard 2004). A similar absor ption pattern was ob - 
served in the 4.8 GHz H2CO transition by Downes et al.1 (Il980l) . 
who associated G23.01-0.41 to a nearby diffuse H 11 region, 
identified through the HI 10a recombination line emission. 

On angular scales of a few arcsec , NH3 (3, 3) VLA 
and CH3CN (6 - 5) NMA observations ( ICodella et alj [19971 



iFuruva et al.ll2008l) revealed the presence of an HMC, with an 
elongated shape and a velocity gradient oriented in the NW- 
SE direction. From the 3 mm dust continuum emission, the 
HMC mass was estimated to be of about 70 Mq (IFuruva et al.l 
2008). These observations suggest that the HMC may be a 
flattened structure rotating about the axis of the '^CO (1 - 
0) bipolar outflow detected at larger scales. The VLA posi- 
tions o f the H2O, OH (Fors ter & Caswefl 1999), and 4.8 GHz 
H2CO dArava et al.ll200ir masers, as wel l as the positions of 
the 6.7 GHz (J. L. CasweU unpu blished; ICaswefl etal] l2000h 
and 12.2 GHz dBrunthaler et al.l [2009[) CH3OH masers mea- 
sured with the Australia Telescope Compact Array (ATCA) and 
Very Long Baseline Array (VLBA), respectively, agree well with 
the position of the HMC. From the large collection of multi- 
wavelength data reported here, the G23. 01-0.41 region can be 
convincingly depicted as an active site of massive star forma- 
tion. 



3. Observations and Calibration 

3.1. VLA observations: 1.3 cm & 3.6 cm continuum 

The source was observed with the VLA[3 at X and K bands in 
both the C-array (in April 2008) and A-array configurations (in 
October 2008). At 3.6 cm in both configurations and at 1.3 cm 
in the VLA-C, the continuum mode of the correlator was used, 
resulting in an effective bandwidth of 172 MHz. At 1.3 cm in 
the VLA-A we used mode "4" of the correlator, with a pair of 
3.125 MHz bandwidths (64 channels) centered on the strongest 
H2O maser line and a pair of 25 MHz bandwidths (8 channels) 
sufficiently offset from the maser lines to obtain a measurement 
of the continuum emission. The two bandwidths centered at the 
same frequency (measuring the two circular polarizations) were 
averaged. 

At X band, 3C286 (5.2 Jy) and 3C48 (3.1 Jy) were used as 
primary flux calibrators, while 1832-105 (1.4 Jy) was the phase 
calibrator For the K-band observations, the primary flux calibra- 
tor was 3C286 (2.5 Jy), the phase calibrator 1832-105 (1.0 Jy), 
and the bandpass calibrator (for the VLA-A data only) 1733-130 
(3.7 Jy). 

The data were calibrated with the NRAO AIPS software 
package using standard procedures. Only for the VLA-A data 
at 1.3 cm, several cycles of self-calibration were applied to the 
strongest maser channel, and the resulting phase and amplitude 
corrections were eventually transferred to all the other line chan- 
nels and to the K-band continuum data. This procedure resulted 
in a significant (at least a factor 2) improvement of the signal-to- 
noise ratio (SNR). 

Natural-weighted maps were made with task IMAGR of 
AIPS both for the continuum and the line data. Simultaneous 
observation of the line and continuum emission at K band made 
it possible to obtain the relative position of the continuum im- 
age with respect to the H2O maser spots with great precision 
(~ O'.'Ol). Finally, we identified the maser spots that appear in 
both our VLBA and VLA images and, using the accurate abso- 
lute position information obtained from the VLBA observations, 
we re-centered the VLA spots and continuum image. The correc- 
tion thus applied turns out to be < 20 mas in each coordinate. In 
conclusion, we believe that the absolute astrometrical precision 
of the continuum maps presented in this paper is of ~ O'.'Ol. 



' The reported values of mass and luminosity have been cor- 
rected taking into acco unt the accurate distance recently measured by 



iBrunthaler etaP ( l2009l) . 



^ The VLA is operated by the National Radio Astronomy 
Observatory (NRAO). The NRAO is a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 



A. Sanna et al.: VLBI study of maser kinematics in high-mass SFRs. II. G23. 01-0.41 



3 



3.2. Maser VLBI 

We conducted VLBI observations of the H2O and CH3OH 
masers (at several epochs), and of the OH maser (at a single 
epoch) toward G23.01-0.41 in the K, C, and L bands, respec- 
tively. In order to derive the maser absolute position, we used 
phase-referencing by fast switching between the maser target 
and the calibrator JI825-0737. This calibrator has an angu- 
lar offset from the maser source of 2°.6 and belongs to the list 
of sources defining the International Celestial Reference Frame 
(ICRF). Its absolute position is known to better than +2 mas 
and its flux measured with the VLB A at S and X bands i s 
68 and 151 mJy beam respectively dFomalont et alj l2003h . 
Five fringe finders (J1642+3948; J1751-I-0939; J1800-I-3848; 
J2101+0341; J2253-i-1608) were observed for bandpass, single- 
band delay, and instrumental phase-offset calibration. Data were 
reduced with AIPS following the VLBI spectral line procedures. 
Paper I describes the general data calibration procedures. 

3.2.1 . VLBA observations: 22.2 GHz H2O masers 

We observed the HMSFR G23.01-0.41 (ti-acking center: 
7;.A.(J2000) = 18''34'"4a:39 and Dec.(J2000) = -09°00'38'.'5) 
with the VLBA0 in the 616 - 523 H2O transition (rest frequency 
22.235079 GHz). The observations (program code: BM244) 
consisted of 4 epochs: April 17, June 29, and September 22, 
2006, and January 17, 2007. During a run of about 6 h per epoch, 
we recorded the dual circular polarization through a 16 MHz 
bandwidth centered on a LSR velocity (Vlsr) of 77.0 km s ' . 
The data were processed with the VLBA FX correlator in 
Socorro (New Mexico) using an averaging time of 1 s and 1024 
spectral channels. The total-power spectrum of the 22.2 GHz 
masers toward G23.01-0.41 is shown in Fig.[T](top panel). This 
profile was obtained by averaging the total-power spectra of all 
VLBA antennas, weighting each spectrum with the antenna sys- 
tem temperature (T^,^). 

The natural CLEAN beam was an elliptical Gaussian with 
a FWHM size of about 1.3 mas x 0.4 mas at a P. A. of -19° 
(east of north), with little variations from epoch to epoch. The 
interferometer instantaneous field of view was limited to about 
2'.'7. At each observing epoch, using an on-source integration 
time of about 2.5 h, the effective rms noise level of the channel 
maps (cr) varied in the range 0.006-0.03 Jy beam"' . The spectral 
resolution was 0.2 km s"' . 

3.2.2. VLBA observations: 1.665 GHz OH masers 

We observed the HMSFR G23.01-0.41 (tracking center: 
R.A.{i2QQQ) = 18''34'"40^:27 and Dec.(J2000) = -09°00'37'.'5) 
with the VLBA in the ^113/2 J = 3/2 OH transition (rest 
frequency 1.665401 GHz), on April 27, 2007 (program code: 
BM244N). During a run of about 6 h, we recorded the dual circu- 
lar polarization through two bandwidths of 1 MHz and 4 MHz, 
both centered at a LSR velocity of 70.0 km s ' . The 4 MHz 
bandwidth was used to increase the SNR of the weak L-band 
signal of the continuum calibrator. The data were processed with 
the VLBA FX correlator in two correlation passes using 1024 
and 512 spectral channels for the 1 MHz and 4 MHz bands, re- 
spectively. In each correlator pass, the data averaging time was 
2 s. The Tjyj-weighted mean of antenna total-power spectra for 
the right and left circular polarizations at 1.665 GHz are shown 
in Fig. [T] (bottom panel). After removing a first order baseline, 



3 The VLBA is operated by the NRAO. 



the 1 .665 GHz OH spectra of both circular polarizations present 
similar characteristics, i.e. several maser emission components 
superposed on a broad absorption feature (see Sect. 14.3b . 

The natural CLEAN beam was an elliptical Gaussian with a 
FWHM size of 19 mas x 10 mas at a PA. of 3°. The interfer- 
ometer instantaneous field of view was limited to about 18'.'5. 
With an on-source integration time of about 1.9 h, the effective 
rms noise level of the channel maps was about 0.02 Jy beam"'. 
The 1 MHz band spectral resolution was 0.2 km s"' . The visibil- 
ity phase of both circular polarizations was calibrated using as 
phase-reference the brightest maser channel of the right circular 
polarization. 

3.2.3. EVN observations: 6.7 GHz CH3OH masers 

We observed the HMSFR G23.01-0.41 (tracking center: 
R.A.ii2Q0Q) = 18''34"W.39 and Dec.(J2000) = -09°00'38'.'5) 
with the European VLBI Network (EVNf| in the 5i - 60A+ 
CH3OH transition (rest frequency 6.668519 GHz). This work 
is based on 3 epochs (program codes: EM061, EM069), sep- 
arated by about 1 yr, observed on February 27, 2006, on 
March 17, 2007, and on March 16, 2008. At the first two 
epochs, antennas involved in the observations were Cambridge, 
Jodrell2, Effelsberg, Hartebeesthoek, Medicina, Noto, Torun, 
and Westerbork. Since the longest baselines involving the 
Hartebeesthoek antenna (e.g., Ef-Hh baseline about 8042 km) 
heavily resolve the maser emission and do not produce fringe- 
fit solutions, the Hartebeesthoek antenna was replaced with the 
Onsala antenna in the third epoch. During a run of about 6 h per 
epoch, we recorded the dual circular polarization through two 
bandwidths of 2 MHz and 16 MHz, both centered on a LSR ve- 
locity of 77.0 km s"' . The 16 MHz bandwidth was useful to 
increase the SNR of the weak continuum calibrator. The data 
were processed with the MKIV correlator at the Joint Institute 
for VLBI in Europe (JIVE - Dwingeloo, The Netherlands) us- 
ing an averaging time of 1 s and 1024 spectral channels for both 
observing bandwidths. The Effelsberg total-power spectrum at 
6.7 GHz toward G23. 01-0.41 is shown in Fig. [T] (middle panel). 

The natural CLEAN beam was an elliptical Gaussian with a 
FWHM size of about 13 mas x 5 mas at a PA. of 30°, slightly 
varying from epoch to epoch. The interferometer instantaneous 
field of view was limited to about 9'.'2. At each observing epoch, 
using an on-source integration time of about 2.2 h, the effec- 
tive rms noise level of the channel maps varied in the range 
0.008-0.3 Jy beam"'. The 2 MHz band spectral resolution was 
0.09 km s"' . 



4. Results 

The present section reports our main results for the radio contin- 
uum and individual maser species in G23.01-0.41. We describe 
in Paper I the criteria used to identify, derive parameters (posi- 
tion, intensity, flux and size), and measure proper motions of in- 
dividual masing-clouds. In the following, the term "spot" refers 
to maser emission on a single velocity channel, whereas the term 
"feature" indicates a collection of spots emitting at similar posi- 
tion over contiguous channels (i.e. an individual masing-cloud). 



The European VLBI Network is a joint facility of European, 
Chinese and other radio astronomy institutes founded by their national 
research councils. 
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4.1. 1 .3 cm & 3.6 cm continuum 

Our measurements of the 1.3 and 3.6 cm continuum emission 
in both VLA-A and VLA-C configurations are summarized in 
Table [T] While the most extended VLA configuration partly re- 
solved the radio continuum at 1.3 cm (missing about half of the 
flux measured in the compact configuration), the fluxes at 3.6 cm 
of both configurations are comparable. Taking into account the 
shape of the VLA-A beam (HPBW of a.'13 x C'DQ at a P. A. 
of -13°), the 1.3 cm continuum emission appears slightly elon- 
gated in the E-W direction (see Fig.O. Fig.|2]shows that a spec- 
tral index (a) of 0.9 is derived from our highest-flux measure- 
ments. 

Previous observations at 3 and 6 cm were performed us- 
ing the ATCA and the VLA-BnA configuration, respectively 
(iForster & Caswellll2000t lArava et al.ll2008l) . iForster & Caswelll 
OOOtf) reported a flux at 3 cm of 0.3 7±0.12 mJy with an HPBW 
of I'.'S X r.'5, an d iArava etalj(l2008h derived a 5cr upper limit of 
2.4 mJy in the 6 cm band (HPBW of 1'.'6 x I'.'l at a PA. of -7°). 

4.2. 22.2 GHz H2O masers 

We imaged the wh ole range of LSR velocities, from 68 to 
87 km s ' , where IForster & Caswe detected water 

maser emission. Using the VLA-C they identified two regions 
of maser activities, each elongated in the north-south direc- 
tion by about 1", and separated by about 40". Our field of 
view covers the region associated with their brightest spots 
and the absolute position we have determined for the reference 
maser sp ot is offset by about I'.'S from the reference position of 
IForster & C aswell (1999). 

We have identified 86 distinct water maser features dis- 
tributed within an area of about 1" x O'.'S. They are mainly 
grouped in three clusters to the NE, SW, and towards the cen- 
ter of the plotted area, which are identified in Fig.|3t with labels 
"A", "B", and "C", respectively. The individual features proper- 
ties are presented in Table|2] Their intensities range from 0.07 to 
202.89 Jy beam"'. Only 28 features (33% of the total) persisted 
over at least 3 epochs, 15 of which lasting 4 epochs. The spread 
in LSR velocities ranges from 89.4 km s ' , for the most red- 
shifted feature (# 57), to 57.0 km s"' , for the most blueshifted 
one (# 22). The clusters "A" and "B" are observed at a similar an- 
gular separation of about 0'.'5 (i.e. 2 x lO"* AU) from the cluster 
"C" and both have a threadlike N-S morphology extending up 
to 1 17 mas. Water maser emission from cluster "A" (15 features) 
and "B" (6 features) presents different properties: the cluster "B" 
consists of faint (weaker than 0.49 Jy beam '), short-living (one 
epoch) features, whereas the cluster "A" consists of bright (up 
to 22.21 Jy beam"') and persistent emission (67% of the fea- 
tures lasted at least 3 epochs and 8 features persisted 4 epochs). 
The cluster "C" has an arc-like morphology, of about 200 mas 
in size, and extends on top of the 1.3 cm continuum emission. 
Water maser emission from this region (61 features) is the most 
variable: only 25% of the observed features persisted for at least 
3 epochs and only 4 features lasted 4 epochs. This region has 
also undergone a powerful maser flare (feature # 1 in Table |2]l 
observed during the last epoch (Fig. [1] upper panel) with an in- 
crease in intensity of about 200 times. The water masers abso- 
lute distribution and line of sight (l.o.s.) velocities are plotted in 
Fig. [3^. Using a time baseline of 9 months, we measured relative 
transverse velocities of individual features with a mean accuracy 
of about 30%. The magnitude of relative proper motions ranges 
from 6.0+3.7 km s"' , for feature # 21, to 59.3±4.7 km s"' , for 
feature # 6. Positions and velocities are relative to the isolated 



and compact feature # 24. Our measurements of absolute proper 
motions have been corrected by the velocity vector i - -37 
and jjy - -90 km s"' (where the X and Y indices refer to the 
motion components toward the east and the north, respectively) 
representing t he Galactic peculia r motion of the HMSFR as mea- 
sured by Brun thaler et al.l(l2009l) . Mean uncertainties in the mag- 
nitude of the absolute motions of water maser features are about 
40%. The relative and absolute proper motions of water masers 
are plotted in Figs.jS^ andjSj^, respectively. Absolute velocities 
of maser features belonging to cluster "C" show a regular vari- 
ation of orientation with position, rotating from north to south 
through the west across the arc-like distribution of features, and 
it appears that they are diverging from a point to the east of the 
maser cluster 

4.3. 1.665 GHz OH masers & absorption 

We imaged the whole range of LSR velocities, from 65 to 
81 km s~', where 1.6 65 GHz OH maser emission was detected 
by IForster & Caswelll ([l999). Using the VLA A-B hybrid con- 
figuration they found OH maser emission scattered over a re- 
gion of about 2" x 1". The absolute position of our reference 
maser spot is offset by about 0'.'9 from the reference position of 
Forster& Caswell (1999). 

Although the whole distribution of hydroxyl masers is scat- 
tered across a region of about 0'.'8 x 0'.'4, half of maser features 
are clustered into a small area, 0'.'06 x 0'.'12, eastward of the 
1.3 cm continuum peak (Fig. [3};). We have identified 24 dis- 
tinct OH maser features and individual features properties are 
presented in Table [3] Positions are relative to the isolated and 
bright feature # 1. The absolute distribution and l.o.s. velocities 
of hydroxyl masers are plotted in Fig. |3};. The LSR velocities 
range from 79.7 km s"' , for the most redshifted feature (# 13), 
to 60.5 km s"' , for the most blueshifted one (# 23). The mean 
FWHM line width of individual maser features is 0.3 km s ' . 
The l.o.s. velocities are blueshifted toward the radio continuum 
and closer to V,,., at larger distances. An ordered velocity shift 
from 73.0 to 62.6 km s ' (features with label number 8, 11, 
15, 12, 17, 20, 16, 22) is observed in an elongated structure ex- 
tending N-S for 33 mas 150 AU) over the central region. 
Maser intensities range from 0.12 to 2.39 Jy beam"' with the 
brightest features showing blueshifted l.o.s. velocities. We have 
identified 5 Zeeman pairs satisfying the condition that the circu- 
lar polarization components coincide positionally within 3cr and 
are separated in velocity by more than their typical line width of 
0.3 km s"' , corresponding t o magnetic fiel ds less than 0.5 mG 
(e.g., inset in Fig. [3};; see also lFish & Reidl2006i) . The l.o.s. mag- 
netic field strength ranges from -5.8 to -1-3.6 mG. OH maser 
total-power spectra show an absorption feature superposed to 
the maser emission (Fig. [1] bottom panel). The absorption fea- 
ture extends over about 40 km s~' , from about 50 to 90 km s"' , 
with two deeps at the LSR velocities of 58.5 and 77.3 km s"' . 
Absorption is not visible on the cross-power spectra of the short- 
est VLBA baseline (KP-PT ^ 400 km). Assuming the absorp- 
tion is resolved out, this allows to set a lower limit of about 0'.'09 
to the size of the absorption region. 

4.4. 6. 7 GHz CH3OH masers 

We imaged the whole range of LSR velocities, from 69 to 
84 km s"' , where 6.7 GHz mas er emission was prev iously de- 
tected towards G23.01-0.41 by ICasweU et alj ( Il995h using the 
Parkes 64-m telescope. The ATCA position of the 6.7 GHz 
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maser peak reported bv lCaswell et alJ (l2000l) falls inside the area 
of our 6.7 GHz maser detections. 

Methanol maser emission at 6.7 GHz is distributed within a 
region of about 0'.'6 x 0'.'6 about the 1.3 cm continuum emission 
(Fig.|3j)). We have identified 81 distinct maser features that are 
grouped depending on the l.o.s. velocity: redshifted features lie 
to the south of the radio continuum and blueshifted ones to the 
north. The individual feature properties are presented in Table|4] 
Maser LSR velocities range from 83.2 km s"' , for the most red- 
shifted feature (# 38), to 67.8 km s"' , for the most blueshifted 
one (#81). The average offsets in l.o.s. velocities (with respect to 
Vsys) of the maser features belonging to the southern and north- 
ern regions of the radio continuum are +4A and -2.1 km s"', 
respectively. Maser intensities are on average very bright (30% 
of the features have peak emission greater than 10 Jy beam"'), 
but range from 0.10 to 146.03 Jy beam"'. 67 features (83% of 
the total) persisted over the 2 yr interval covered by our 3 ob- 
serving epochs. For each maser feature. Table |4]reports the mean 
brightness variability defined as the ratio between the variation 
of the brightness (of the strongest spot) and the average of the 
minimum and maximum brightness. The mean brightness vari- 
ability of methanol maser features was about 45% (compared to 
only 20% for the source G16.59-0.05, Paper I), varying from 
a minimum of 5%, for feature # 80, to a maximum of about 
136%, for feature # 48. These estimates of brightness variability 
should be taken as upper limits since they include effects from 
the slightly varying beam shape among different epochs as well 
as amplitude calibration uncertainties. Positions are relative to 
the structurally-stable, compact, and bright feature #10, whose 
intensity varies less than about 6%. Methanol maser positions 
and l.o.s. velocities are plotted in Fig.[3]3. We have calculated the 
geometric center of all the detected features (hereafter "barycen- 
ter", identified with label S in Table |4]l, whose position is indi- 
cated by a star in Fig.[3j5. In deriving the "barycenter" position, to 
reduce the weight of clustered maser features, uniform weight- 
ing has been used, weighting feature positions with the recipro- 
cal of the local feature density. Note that, since only a subset of 
6.7 GHz maser features had a stable (spatial and spectral) struc- 
ture over time, the "barycenter" is not a good point to refer the 
internal motions. We have also calculated the geometric center 
(hereafter "center of motion", identified with label # in TableHJ 
of features with a stable spatial and spectral structure persisting 
over the 3 observing epochs, and refer our measurement of trans- 
verse motions to this point (Fig.|6^). Using a 2 yr time baseline, 
the relative transverse velocities of individual maser features are 
determined with a mean accuracy of about 30%. Proper motions, 
derived only for maser features with a stable structure, are cal- 
culated by performing a linear fit of feature position vs. time. To 
establish the extent of accelerated motions, an analysis of devia- 
tions from linear motion is postponed after an upcoming fourth 
epoch. The magnitude of relative proper motions ranges from 
1.7+1 km s"' , for feature # 17, to 16.0+1.6 km s"' , for fea- 
ture #15. Transverse velocities of maser features have a mean 
value of 7.2 km s"' , about 2-3 times the spread of l.o.s. ve- 
locities. We observe a cluster of maser features connected by a 
bridge of extended weaker emission (110 mas x 60 mas), emit- 
ting over a narrow velocity range of a few km s"' (Fig. ^p). 
Since in this region spatial blending of nearby maser centers 
limits the precision of positions derived with Gaussian fitting, 
proper motions from this area were sufficiently accurate only for 
the brightest, more isolated features (# 1, 4, 9, 17, 18, 19). 

Looking at Fig. [3}', 6.7 GHz maser features seem to trace a 
bipolar, funnel-like structure, centered to the east of the 1.3 cm 
continuum peak. Maser proper motions presented in Fig. [6^ de- 



scribe a smoothly varying velocity field, roughly symmetrical 
with respect to the barycenter of the spatial distribution. The 
regularity and near symmetry of the maser velocity distribution 
makes us confident that the calculated proper motions are repre- 
sentative of the true motions of the methanol masing-clouds. 

5. Discussion 

5.1. Nature of the radio continuum 

The continuum source appears slightly resolved at the VLA-A 
1.3 cm angular resolution of ~ O'.'l (see Figs.|4H6]l. The fluxes 
of the continuum source at 1.3 and 3.6 cm (see Table [TJ are con- 
sistent with free-free emission from an UCH ii region with ra- 
dius of about 450 AU ionized by a Lyman continuu m of 1-3 x 
lO"**" s"', corresponding to a ZAMS Bl star ( Schram l & Mezge^ 
1 19691: iPanagial 11973). The spectral index a = 0.9 between 
3.6 and 1.3 cm could be interpreted with the emission becom- 
ing optically thin at /i $ 1 cm. 

An alternative interpretation for the weak continuum emis- 
sion is in terms of a shock in duced ion ization in a thermal 
jet. Following the calculation bv lAngladal (Il996l) . for optically 
thin emission, one has: Fd^ - 10^'^ (£l/4n)P, where F is the 
measured continuum flux in mJy, P is the outflow momentum 
rate in M0 yr"' km s"' , Q is the jet solid angle in sr, and d 
is the source distance in kpc. Using the flux of 1.4 mJy mea- 
sured at 1.3 cm with the VLA-C, and the accura te distance of 
4.6 kpc to G23.01-0.41 bv lBrunthaler et al.l (l2009h . one derives: 
P - 10"^ (0/47r)"' M© yr"' km s~' . The momentum rate thus 
depends on the estimate of the jet collimation factor (Q/47r)"'. 
If the jet structure is not resolved but only a bright knot along 
the jet axis is observed, the jet solid angle can be estimated from 
the angular size of the knot if the position of the YSO power- 
ing the jet is known. In the next section we discuss the maser 
information with the aim to constrain the YSO position. 



5.2. Water masers tracing a wind/jet system 

Current models explain H2O maser excitation by coUisional 
pumping with H2 molecules within hot (> 400 K) shocked lay- 
ers of gas b ehind both high-ve locity (> 50 km s"') dissocia- 
tive (J-tvpe: lElitzur et all 19891) and slow (< 4 5 km s"') non- 
dissociative (C-tvpe: lKaufman & Neu feld* 1996) shocks, propa- 
gating in dense regions (H2 pre-shock density in the range 10^- 
10^ cm"^). Following these models, the cluster "C" of water 
masers witnesses the presence of dense shocked gas near the 
continuum source. The higher variability of the bright water 
masers from this cluster, compared with clusters "A" and "B", 
could indicate that the cluster "C" pinpoints the most active site 
of the region. The arc-like distribution of maser features in this 
region, together with their fast (as high as 55 km s"' ) and diverg- 
ing proper motions (see Fig. [Sji), lead us to think that this water 
maser emission traces a shock front expanding from a center lo- 
cated to the east of the 1.3 cm continuum peak. The barycenter 
of the funnel-like structure of the 6.7 GHz masers (the star in 
Fig. |3]l falls also in this area, and we speculate that is the place 
of the massive YSO(s), possibly a multiple system, exciting and 
driving both maser emissions. In the following discussion, we 
will use the 6.7 GHz maser barycenter as the best guess for the 
YSO(s) position. 

If water maser emission is excited in a jet from a YSO, 
knowing the average distance of water masers from the YSO 
and the average maser velocity, one can estimate the momen- 
tum rate of the jet assuming that is momentum driven. Using 
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the model-predicted H2 pre-shock density of 10^ cm ^, the mo- 
mentum rate in the water maser jet is given by the expression: 
P = 1.5 X 10-^YIqRI^^(^1/47t) Mgyr-i kms"', where Vjo 
is the average maser velocity in units of 10 km s ' , Rioo is 
the average distance of water masers from the YSO in units of 
100 AU, and O is the solid angle of the jet. This expression 
has been calculated by multiplying the momentum rate per unit 
surface transferred to the ambient gas (n/^, m^T V^), by a factor 
QR^, assuming that the jet is emitted from a source at a distance 
R from the masers within a beaming angle Q. 

We hypothesize that a wind emitted from the YSO impacts 
against the dense circumstellar gas, and excites and drives the 
motion of water masers of cluster "C". Water masers of cluster 
"C" have an average (sky-projected) distance from the YSO of 
» 200 AU, and their average absolute velocity is 40 km s"' . 
Assuming an isotropic wind (Q = An), and considering that the 
derived sky-projected distance should be taken as a lower limit 
to the true distance, we derive a value for the momentum rate 
P > 0.1 Mq yr"' km s"' . This value is one order of magnitude 
higher than the typical momentum rate of the stellar wind from 
a ZAMS O star of 10"^ Mq yr"' km s ' (Markova et al. 20051), 
and might be explained if the massive YSO was still accreting 
mass and a fraction of mass approaching the YSO was deflected 
into t he wind, as predicted by X-wind (Shu et al. 2000) and disk- 
wind (Konigl & Pudritz"2000) models. 

Assuming that the position of the YSO is close to the geo- 
metric center of the (methanol and water) maser distributions, 
the 1.3 cm continuum source is found off'set from the YSO by 
about 60 mas to the SW. Such an offset is significantly larger 
than the accuracy (w 10 mas) of the absolute position of the 
continuum source, which leads us to think that the best interpre- 
tation is in terms of a thermal jet rather than an UCH 11 region. 
We describe a case similar to what observed for the well-studied 
source IRAS 20126+4104, where the YSO position is clearly 
off'set of the peaks of the VLA 3.6 cm continuu m emission, 
tracin g shock induced ionization in a thermal jet (.Hofner et al.l 
I2007h . Towards G23.01-0.41, the VLA-A slightly resolves the 
1.3 cm emission, which indicates that the source size is less than 
the VLA-A beam at 1.3 cm of w 110 mas. Then, at the dis- 
tance of the 1 .3 cm continuum peak from the YSO of x 60 mas, 
the jet has to be collimated within a cone of semi-aperture 
6 < 0.5 X 110/60 = 0.9 rad, corresponding to a jet solid angle 
f2 < 2.5 sr. Using this value in the expression for the momentum 
rate derived in Sect. 15. II one has P > 5 x 10"^ M0 yr ' km s"' , 
which is consistent with the value of the momentum rate derived 
above for the YSO wind. 

Figure |3^ shows that water masers not belonging to cluster 
"C" are distributed to the NE (mainly in cluster "A") and to the 
SW (mainly in cluster "B") of the putative location of the YSO. 
They could trace shocked gas along the path of the bipolar jet, 
emerging from the YSO and interacting with the surrounding 
molecular core. Such a jet could be the engine of the massive 
'-CO (1-0) bipolar outflow, elongated in the NE-SW direc- 
tion, observed on arcsec scales (see the sketch in Fig. [8]l. This 
interpretation explains the absolute velocities of water masers 
in cluster "A" oriented towards NE (see Fig. The same jet 
should be responsible for the excitation of the continuum source 
and all the water maser features separated from the YSO by more 
than w 0'.'2. To intercept all the water masers found at larger dis- 
tance from the YSO, the jet semi-aperture has to be larger than 
> 0.25 rad, which is consistent with the condition 9 < 0.9 rad 
from the upper limit to the angular size of the continuum source. 

Water masers in cluster "A" move with an average velocity 
of 20 km s ' and are separated from the YSO by » 2300 AU. 



Using these values in the expression for the momentum rate in 
the maser jet, one finds: P - 3.2 (Q/47r) M0 yr"' km s"' . If the 
same jet is responsible for exciting the continuum source and 
driving the water masers of cluster "A", the jet solid angle (Q.) 
can be determined by requiring that the momentum rate in the 
maser jet equals that from the continuum emission: 



3.2(Q/47r) 



10"^(Q/47r)"' 



From this equation one derives Q. = 0.7 sr corresponding to a 
semi-opening angle of the jet - 0.5 rad. This value, derived 
from measurements of the continuum flux and the maser kine- 
matics of cluster "A", is consistent with the independent condi- 
tion 0.25 < 6 < 0.9 owing to the spatial distribution of all 
detected water masers and to the size of the continuum source. 
That supports our result. From the derived value of the jet solid 
angle, the jet momentum rate is P = 0.2 M0 yr"' km s"' . 

This value for the jet momentum rate agrees with that 
(0.1 Mq yr"' km s"' ) of the dominant (SW redsh ifted) lobe of 
the lar ge-scale '^CO (1-0) outflow reported bv iFuruva et al.l 
(l2008h . These evidence support the interpretation that the NE- 
SW distribution of water masers traces the jet at the base of the 
massive molecular outflow observed at arcsec scales. Note also 
that the momentum rate in the jet is comparable with the mo- 
mentum rate (0.1 Mq yr"' km s"' ) derived for the wide-angle 
wind driving the water masers of cluster "C", which presents 
the case that the jet itself is powered by the YSO wind. At a 
distance of w 300 AU from the YSO, i.e. at the position of the 
1.3 cm source, most of the momentum of the wind would be 
efficiently collimated within the jet, with a collimation factor 
(Q/47r)"' w 20. For comparison, values of collimation factors 
of the order of 10 are typical for well-studie d radio jets in lo w- 
and intermediate-mass star-forming regions (lAngladalll996h . If 
the high momentum rate of the wind could hint to the YSO be- 
ing still accreting mass, the efficient collimation of the wind into 
a jet could indicate that the accretion onto the YSO is mediated 
through an accretion disk. 

5.3. Methanol maser environment 



In Sect. 15.21 we have used the symmetries in the spatial distribu- 
tions of the 6.7 GHz methanol and 22 GHz water masers to con- 
strain the YSO position, and postulated that the YSO is located 
close to the geometric center of both maser species distributions. 
Figure m shows that 6.7 GHz and 22.2 GHz masers emerge from 
nearby but different positions around the putative location of 
the massive YSO. Figures ^ and [6^ show that water masers 
move significantly faster than methanol masers, and the average 
direction of motion of the two maser species is also diff'erent. 
Therefore, the comparison of the spatial and velocity distribu- 
tion suggests that the two maser species, although associated 
with the same massive YSO, originate from distinct environ- 
ments, characterized by a diff'erent kinematics and, likely, as pre- 
dicted by water and methano l maser excitation models, by differ- 
ent ph ysical conditions (e.g.. lMoscadelli et alj2007t Goddi et al] 
l2007h . CH3OH excitation models (ICragg et al.1 l2005h predict 
that Class II methanol maser emission is produced by radia- 
tive pumping in cool (a: 30 K), or moderately warm regions 
(< 200 K), with H2 densities in the range 10''-10'^ cm"^. Even 
though comparably high densities are required for both water 
and Class II methanol maser action, strong 6.7 GHz methanol 
masers can be produced at gas temperature significantly lower 
than that of the shocked layers of gas emitting water masers. 

Looking at Fig.|6^, one notes that the distribution of 6.7 GHz 
maser velocities is not isotropic about the maser barycenter 
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Taking the maser barycenter as the origin of a system of polar co- 
ordinates, (sky-projected) maser velocities with a dominant ra- 
dial component appear to concentrate at a P. A. different from that 
where maser velocities are predominantly azimuthal. Figure [7] 
shows the distribution of the angle between the maser proper 
motion and the corresponding position vectors, as a function of 
both the radial distance and the P. A. of the maser feature. Most of 
maser velocities form an angle with the radius (outward from the 
maser barycenter) varying in the range from -20° to 120° (evalu- 
ated positive counterclockwise for the observer). While the dis- 
tribution of the velocity-position angle appears to be rather uni- 
form with the maser radial distance, all the azimuthal velocities 
(forming with the radius an angle within 90° + 30°) concentrate 
in a well defined PA. range, 20° < P.A. < 70°. 

The observed distribution of 6.7 GHz maser velocities can 
be interpreted in terms of a combination of expanding and rotat- 
ing motions. Let us consider a rotating structure about an axis 
inclined with respect to the l.o.s., and indicate with and y^x 
the projections on the plane of the sky of the rotation axis and the 
line perpendicular to that, respectively (see the sketch in Fig.|8]l. 
Points of the rotating structure belonging to the plane contain- 
ing the rotation axis and the l.o.s. lie on Xax and have rotation 
velocities parallel to y^x. Therefore, it is possible to identify the 
projection onto the plane of the sky of the rotation axis looking 
for the direction in the sky to which the proper motions are per- 
pendicular. The upper plot of Fig.|7]suggests that the azimuthal 
6.7 GHz maser velocities can result from rotation around an axis 
whose Xax crosses the maser barycenter and is oriented at P.A. 
within 45° + 25°. Note that this Xax crosses also the 1.3 cm con- 
tinuum peak (at a P.A. of 225°) and agrees with the orientation 
of the 22 GHz water maser jet discussed in Sect. 15.21 Then, there 
are indications that 6.7 GHz masers can emerge from a structure 
rotating around the jet/outflow axis. 

Many of the observed 6.7 GHz maser velocities show a com- 
ponent directed radially outward from the maser barycenter, sug- 
gesting that the gas traced by 6.7 GHz masers is also expanding 
away from the YSO. Looking at Fig. [6^ and the upper plot of 
Fig. |7] one notes that a group of approximately radial velocities 
(directed within 10° from the local radius) concentrates at P.A. 
between 150° - 160° and 330° - 340°, that is along a direction 
approximately perpendicular to Xax. If the rotation axis is close 
to the plane of the sky, the rotation velocities of the spots be- 
longing to the plane containing the rotation and yax axes should 
be almost parallel to the l.o.s.. Then, the sky-projected velocities 
along the yax axis are expected to show small azimuthal compo- 
nents and to be mainly radial if the gas is also flowing away from 
the YSO, in agreement with what we actually observe. 

The same assumptions that the 6.7 GHz masers trace gas 
flowing away from a center and rotating about an axis at small 
inclination with the plane of the sky, allows also to easily inter- 
pret the observed distribution of l.o.s. velocities. Material on the 
near and far side of the maser structure would move towards and 
away from the observer, respectively. Along the sky-projected 
rotation axis (i.e. Xax), the observed blue- and redshifted maser 
velocities towards the NE and S W, respectively, can be explained 
if we admit that maser emission to the NE and SW originates 
from material on the near (approaching) and far (receding) side 
of the structure, respectively. Along the direction perpendicular 
to the sky-projected rotation axis (i.e. yax), the distribution of 
l.o.s. velocities would reflect the sense of rotation, with blue- 
and redshifted maser velocities to the NW and SE, respectively, 
in agreement with the observations (Fig.[8]l. 

If the rotation axis is slightly inclined with respect to the 
plane of the sky, the small group of maser features along the yax 



axis with approximately radial velocities (negligible azimuthal 
components), should be found close to the sky. Then, their dis- 
tribution of l.o.s. velocities can be used to constrain the central 
mass, in the hypothesis of centrifugal equilibrium. Their aver- 
age distance from the maser barycenter (Fig. |71 lower plot) is 
~ 1000 AU, and the average of the absolute value of their l.o.s. 
velocities is a: 4 km s"' . Using these values, a central mass of 
~ 20 Mq is derived, corresponding to an early B - late O type 
ZAMS star. 

This discussion suggests that the geometrical properties and 
the amplitude of the observed 6.7 GHz maser velocities can 
be explained in terms of rotation and expansion about a mas- 
sive YSO. The idea that the methanol-masing gas can partici- 
pate in an overall expan ding motion has also been proposed by 
iBartkiewicz et al. 1 12009 1) on the basis of single-epoch VLBI ob- 
servations of a sample of 31 6.7 GHz maser sites. These authors 
have noted that most maser regions coincide with a 4.5 //m emis- 
sion excess (as it is also the case for G23.01-0.41), which is 
commonly interpreted as a tracer of expanding, shocked, molec- 
ular gas. In the G23. 01-0.41 region, the derived orientation of 
the rotation axis (at P.A. within 45° + 25°) is consistent with 
the scenario where 6.7 GHz CH3OH masers can be tracing the 
internal portions (within a radius of about 2000 AU from the 
YSO) of the rotating core, observed in the NH3 and CH3CN 
lines on arcsec sca le along the SE-NW direction (ICodeUa et alJ 
1997; F uruva et al.k200 8). The pattern of CH3CN l.o.s. velocities 
observed bv lFuruvaet a l. (2008, Fig. 11), mainly blueshifted to 
the NW and redshifted to the SE, agrees with the distribution of 
l.o.s. velocities of 6.7 GHz CH3OH masers at much smaller an- 
gular scale. We note that the 12.2 GHz maser emission measured 
in the G23.01-0.41 region ( Brunthaler et al. 2009) has both po- 
sition and LSR velocity (Vlsr - 81.5 km s"' ) consistent with 
the 6.7 GHz masers (Fig.|4]l. In order to produce bright 6.7 and 
12.2 GHz masers ( Tb > 10'° K), th e Class II methanol maser 
pumping models by Cragg et al. (2005) prescribe a combination 
of warm dust (Td = 175 and 125 K) and cool gas (T^ = 30 and 
50 K), for a gas density of nnj - 10''-10^ cm"-'. These values 
agree well with the parameters of the ammonia core (Tk = 58 K 
and ne, = 6.9 x lO'' cm"^) derived by Codella et al. (I l997h . and 
would marginally support the idea of an association of methanol 
masers with the NH3 core. 

The measured pattern of 6.7 GHz maser proper motions of- 
fers a simple explanation that the CH3CN LS R velocity gradient 
in the toroid mapped by iFuruv a et al.l (1200 8h presents a not so- 
well defined orientation. Gas close to the equatorial plane would 
indeed undertake a complex motion resulting from a combina- 
tion of rotation and expansion. Taking all the 6.7 GHz maser 
features with measured proper motions, the mean absolute value 
(w 4 km s"' ) of the radial velocity components is basically 
the same as that of the azimuthal velocity components, indicat- 
ing that the masing gas expands and rotates at similar speed. 
Assuming that the expansion of the gas traced by the 6.7 GHz 
masers is momentum driven, the momentum rate to accelerate 
the gas can be estimated with the expression P = MV^/R, 
where M is the total mass of the expanding gas, V the av- 
erage expansion velocity and R the radius of the expanding 
sphere (working in spherical symmetry, for simplicity). Taking 
an upper limit to the gas density of ne, = 10^ cm"^, a velocity 
V = 4 km s"' , and an upper limit to the observed (sky-projected) 
distance of 6.7 GHz masers from the YSO of R = 2000 AU, a 
momentum rate P = 3 x lO"-' Mq yr"' km s"' is derived. This 
value is about two orders of magnitude smaller than the momen- 
tum rate of the YSO wind supposed to be driving the expansion 
of water masers of cluster "C", indicating that a negligible frac- 
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tion of the wind, non-collimated into the jet and blowing across 
the equatorial plane, could be responsible for the expansion of 
the gas nearby the YSO. The YSO wind is only one of the possi- 
ble causes for the observed 6.7 GHz maser outward motion. The 
spatial sampling of the measured maser velocities is not suffi- 
ciently complete for a quantitative analysis to effectively con- 
strain the geometry of the masing gas and determine which por- 
tions of it is expanding. 



5.4. A molecular expanding layer in front of the central 
source? 

Figures [3j; and |4] show that the spatial distribution of the 
1 .665 GHz OH masers differs from that of 6.7 GHz CH3OH and 
22.2 GHz H2O masers. The N-S elongated structure traced by 
the OH masers over the 1.3 cm continuum is slightly offset to 
the east of the water maser cluster "C", and is found near to 
the (putative) YSO location. To explain the strongly blueshifted 
1.0. s. velocities (from -4 to -14 km s"' ) of this cluster of OH 
masers, a possible interpretation is that they are seen in the fore- 
ground of the YSO and the radio continuum emission, and that 
they are originating in a layer of gas moving towards the ob- 
server The same wind, which we have conjectured to be driv- 
ing the expansion of water masers in cluster "C", might be re- 
sponsible for accelerating the gas hosting the OH maser struc- 
ture. The O H and CH3OH Class II maser excitation models by 
ICragg et"ai] (2002, Fig. 4) predict strong 1.665 GHz maser ac- 
tion and inhibition of the 6.7 GHz masers when nn, < 10^ cm"^. 
This suggests that the 1.665 GHz OH masers towards the ra- 
dio continuum are tracing lower densities than those required 
for intense 6.7 GHz CH3OH masers. Thus, the properties of the 
elongated structure observed (only) in the 1.665 GHz OH maser 
emission might be explained if the OH masers trace a gas layer 
with nHj < 10'' cm"^, rapidly expanding away from the YSO 
towards the observer. From Zeeman splitting, the l.o.s. magnetic 
field strength in this OH maser layer is estimated to be about 
1 mG (see Table[3] features #11 and 20). This result is in gen- 
eral agreement with a detailed VLBA study of OH maser proper 
motions and polarization properties in several UCH 11 regions by 
iFish & R eid (2007). These authors found that OH masers trace 
the expansion of ionized gas and also established a correspon- 
dence between the gas density and the magnetic field, deriving a 
value of the order of 1 mG for the l.o.s. magnetic field strength 
if nH, = 10'' cm-^ 

We propose that the OH absorption is spatially unre- 
lated to the OH maser emission and may instead be con- 
nected with the large scale dynamics of the region. In the 
NRAO VLA Sky Survey, the G23.01-0.41 region shows dif- 
fuse 1 .4 GHz radio continuum emission extended on scales of 
about r. Toward G2 3 .01-0 .41, using the Effelsberg radiote- 
lescope iDownes et al.l (I198Q) reported both diffuse HllOff re- 
combination line (4.874 GHz) emission and absorption in the 
4.8 GHz H2CO line at various velocities ranging from 55 to 
103 km s"' . Furthermore, a lower limit of a few arcminutes to 
the size of the absorption region is set by the H2CO VLA-D ob- 
servations by lArava etaH (l2008l) . who resolved the absorption. 
The 1.6 GHz OH absorption feature visible in our VLBA total- 
power spectra is mainly blueshifted and that may suggest that 
the absorbing molecular gas stays in front of the extended radio 
continuum emission and is expanding towards the observer (e.g., 
ISzvmczak & Gerardll2004 . 



6. Summary and Conclusions 

Using the VLBI technique, we observed the HMSFR G23.01- 
0.41 in the three most powerful maser transitions: 22.2 GHz 
H20, 6.7 GHz CH3OH, and 1.665 GHz OH. The source G23.01- 
0.41 was also observed with the VLA, detecting faint (* mJy) 
radio continuum emission with 0'.' 1 resolution towards the cen- 
ter of the HMC. Our main conclusions can be summarized as 
follows: 

1. H2O, CH3OH, and OH maser emissions are distributed 
within X! 2000 AU from the center of a HMC powering a 
massive bipolar outflow. The three maser species, although 
associated with the same YSO, present different, although 
complementary, spatial distributions and kinematical prop- 
erties. 

2. Water masers trace fast (20-50 km s"' ) outflows emitted 
from the YSO. An arc-like structure of water masers super- 
posed on the radio continuum marks a fast shock propagating 
through dense gas, and is probably driven by a YSO wind. 
The 1.3 cm continuum source and the two clusters of wa- 
ter masers, aligned along the NE-SW direction detected at 
larger distance from the YSO, are likely tracing the jet driv- 
ing the massive CO outflow observed at larger scales. 

3. Methanol masers present a N-S oriented, funnel-like spa- 
tial distribution with red- and blueshifted features located to 
the south and the north, respectively. Observing 3 different 
EVN epochs spanning 2 yr, we have measured accurate (rel- 
ative errors < 30%) proper motions of 6.7 GHz methanol 
features. The maser transverse velocities range from a few 
to about 15 km s ' . The pattern of 6.7 GHz maser proper 
motions can be interpreted in term of a composition of ex- 
pansion and rotation around a YSO of about 20 Mq), with 
the rotation axis oriented on the sky at similar RA. to the 
axis of the jet/outflow system traced by the water masers. 
It is then plausible that the CH3OH masers trace the internal 
portions of the toroid, elongated along the SE-NW direction, 
observed in the CH3CN and NH3 lines on arcsec scale. 

4. Hydroxyl masers superposed on the radio continuum are 
probably seen in the foreground and expand outward from 
the central source tracing a lower density environment than 
that harboring the methanol and water masers. 

This study demonstrates that multi-epoch VLBI observa- 
tions of different maser species provide information useful to 
explore the complex phenomena occurring at distances of tens 
up to thousands of AU around massive YSOs. The best scien- 
tific return from VLBI maser observations can be expected when 
such observations will be compared to data of new generation 
(sub)millimeter interferometers (such as ALMA). Clearly, ther- 
mal line observations of both outflow(s) and core(s) with an an- 
gular resolution closer to that of VLBI maser data will help to 
solve ambiguities in the interpretation of kinematic structures, 
thanks to a better sampling of the (proto)stellar environment and 
its physical properties. 
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Fig. 1. Total-power spectra of the H2O, CH3OH and OH masers toward G23. 01-0.41. Upper panel: system-temperature (Tjvj) 
weighted average of the 22.2 GHz total-power spectra of the 9 VLBA antennas observing on 2007 January 17 (BR, FD, KP, LA, 
MK, NL, OV, PT, SC). Inset: zoom of the water maser peak. Middle panel: Effelsberg total-power spectrum of the 6.7 GHz methanol 
maser emission on 2007 March 17. Inset: zoom of the methanol maser peak. Lower panel: Tjvj- weighted average of the 1.665 GHz 
total-power spectra of the 9 VLBA antennas observing on 2007 April 27 (BR, FD, HN, KP, MK, NL, OV, PT, SC). Continuous 
(black) and dashed (red) lines are used to distinguish between right (RCP) and left (LCP) circular polarizations, respectively. The 
dotted line crossing the spectra represents the systemic velocity {V sys) inferred from CH3CN measurements. 
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Fig. 2. Spectral energy distribution of the radio continuum in the HMSFR G23.01-0.41. Dots and errorbars report the values and 
the associated eiTors (Icr) of our measurements at 1.3 and 3.6 cm using both the VLA-A and VLA-C configurations (see TablelTJ. 
Taking at each observing frequency the highest-flux measurement, the derived spectral index, indicated by the dotted line, is o- = 
0.9 . 
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Fig.3. Absolute positions and LSR velocities of the three maser species observed in G23. 01-0.41: a) H2O (triangles), b) CH3OH 
(dots), and c) OH (squares). Different colors are used to indicate the maser LSR velocities, according to the color scale on the 
right-hand side of the plot, with green representing the systemic velocity of the HMC. The VLA 1.3 cm continuum emission is 
plotted with dashed contours. Contour levels range from 30 to 90% of the peak emission (0.72 mJy beam"') at multiples of 10%. 
The restoring beam is shown in the lower left corner of the upper panel. The linear scale of the plots is shown in the lower left 
corner of the middle panel. For water masers, letters A, B, and C indicate the clusters defined in Sect. 14.21 For methanol masers, the 
star marks the position of the barycenter of the methanol maser distribution, as defined in Sect. 14.41 For hydroxyl masers, full and 
empty squares denote right and left circularly polarized features, respectively, and numbers close to (some) maser features indicate 
the inferred local value of the magnetic field strength (in mG). Inset: example of Zeeman pair for the feature #11, reporting the 
Gaussian fits to the spectral profiles of the RCP (black line) and LCP (red line) components. The RCP and LCP measurements are 
indicated with full and empty squares, respectively. 
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Fig. 4. Collection of the sub-arcsec observations toward G23.01-0.41. Left panel: map in the NH3 (3,3) line from ICodella et all 
(Il997h . Contour levels range from 20 to 90% of the peak emission (46 mJy beam ') at multiples of 10%. The restoring beam is 
shown in the lower left corner of the panel. The star marks the peak position of the 3 mm continuum emission as determined by 
iFuruva et al.l (l2008l) . Right panel: enlargement of the region over which we have detected maser emission. Red triangles, blue dots, 
and green squares represent respectively H2O, CH3OH, and OH maser positions from our VLBI measure ments. The empty circle 
marks the position of the brightest 12.2 GHz CH3OH maser feature (at Vlsr = 81.5 km s ' ) measured bv lBrunthaler et al"J(l2009h 
with an uncertainty of a few mas. The empty triangle with the errorbars indicates the position (and the associated uncertainty) 
of the 4.8 GHz H2CO maser feature (at Vlsr = 73.6 km s"' ) derived by Araya et al. (2008). Dashed contours show the 1.3 cm 
VLA continuum emission. Contour levels range from 30 to 90% of the peak emission (0.72 mJy beam"') at multiples of 10%. The 
restoring beam is shown in the lower left corner of the panel. 
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Fig.5. 22.2 GHz H2O maser kinematics toward G23.01-0.41. a) positions (triangles) and transverse velocities of the H2O maser 
features relative to the feature # 24 (indicated with the vertex-connected symbol). Colored cones are used to show both the direction 
and the uncertainty (cone aperture) of the proper motion of maser features. The proper motion amplitude scale is given by the black 
arrow on the bottom right corner of the panel. Proper motions of the maser features # 30, 38, and 62 are too uncertain and, thus, 
not plotted. Different colors are used to indicate the maser LSR velocities, according to the color scale on the right-hand side of 
the panel, with green denoting the systemic velocity of the HMC. h) absolute positions and transverse velocities of the water maser 
features. The potted field of view is the same as in the upper panel and symbols have the same meaning as in the upper panel. Dashed 
contours show the VLA 1.3 cm continuum emission. Contour levels range from 30 to 90% of the peak emission (0.72 mJy beam"') 
at multiples of 10%. The restoring beam is shown in the lower left corner of the panel. 
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Fig. 6. 6.7 GHz CH3OH maser kinematics toward G23. 01-0.41. aj absolute positions (dots) and transverse velocities of the CH3OH 
maser features relative to the center of motion (as defined in Sect. 14.4b of the methanol maser distribution (indicated by the cross). 
Colored cones are used to show both the direction and the uncertainty (cone aperture) of the proper motion of maser features. The 
proper motion amplitude scale is given by the black arrow on the bottom right corner of the panel. Different colors are used to 
indicate the maser LSR velocities, according to the color scale on the right-hand side of the plot, with green denoting the systemic 
velocity of the HMC. Dashed contours show the VLA 1 .3 cm continuum emission. Contour levels range from 30 to 90% of the peak 
emission (0.72 mJy beam"') at multiples of 10%. The restoring beam is shown in the lower left corner of the panel, h) Zoom of the 
crowdest region of 6.7 GHz maser emission. The map is obtained from data of the second EVN epoch (2007 March 17), summing 
the maser emission over the velocity range from 73.40 to 74.63 km s"' . Plotted contour levels are spaced by factor of 2, with the 
lowest level corresponding to the 3cr rms of 3.0 Jy beam '. The restoring beam is given in the lower left corner of the panel. The 
numbers near emission peaks are the maser feature labels given in Table |4] 
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Fig. 7. Upper and lower panels show the distribution of the angle between the proper motion and the position vectors of 6.7 GHz 
masers as a function of the maser P. A. and radial distance, respectively. Approximately azimuthal and radial velocities are indicated 
with empty squares and triangles, respectively (velocity-position angle in the range from 60° to 120° and from -10° to 10°, respec- 
tively), whereas filled dots are used for other values of the velocity -position angle. The 6.7 GHz maser barycenter (see Table |4]l is 
taken to be the origin of the (sky-projected) maser positions. The angle between the velocity and the position vectors of a given 
maser feature is evaluated positive counterclockwise and is equal to the difference between the P.A. of the feature velocity and 
position vectors. Maser features with P.A. in the range 180° - 360° have been folded into the P.A. interval 0° - 180° (by subtracting 
180° to their true P.A.). In the upper plot, horizontal dotted and dashed lines mark the boundary and the center, respectively, of the 
interval over which maser velocities are approximately azimuthal. 
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Fig. 8. Schematic cartoon of the main components in the molecular cloud associated with G23.01-0.41 . The drawing is not to scale. 
The components represented are (i) the outflow/jet structure traced by CO and H2O maser emission (black dots) from pc-scale to 
a few thousand of AU; (ii) the toroid structure seen on a large scale in CH3CN and highly excited NH3 molecules, and its inner 
portion traced by CH3OH maser emission (blue and red dots according to the l.o.s. velocity); (iii) the radio continuum emission 
superposed on the arc-Uke H2O maser structure. The star marks the assumed YSO(s) location. 
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Table 2. Parameters of VLBA 22.2 GHz water maser features. For each identified feature, the label (given in Col. 1) increases with decreas- 
ing brightness. Cols. 2 and 3 report the LSR velocity and brightness of the brightest spot (for the epoch specified by the number in brackets). 
Cols. 4 and 5 give the position offset relative to the feature # 24, toward the East and the North directions, respectively. Cols. 6 and 7 report the 
projected components of the feature proper motion relative to the feature # 24, along the East and North directions, respectively. Col. 8 shows the 
relative error of the proper motion amplitude. The absolute position and proper motion of the reference feature # 24 are reported at the bottom of 
the Tabic. 
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78.26 


1.72 


(4) 


435.51 ±0.10 


77.83 ±0.13 


6.0 ±3.6 


-3.4 ± 4.0 


54% 


26 


65.63 


1.64 


(1) 


-66.53 ±0.14 


-60.98 ±0.11 


-29.5 ± 4.7 


12.9 ±4.1 


14% 


27 


68.36 


1.52 


(1) 


-64.97 ±0.13 


-60.60 ±0.12 








28 


81.21 


1.43 


(3) 


-62.98 ± 0.08 


-77.14 ±0.09 








29 


79.95 


1.41 


(4) 


-63.87 ±0.13 


-163.86 ±0.10 


-18.6 ±7.3 


-35.6 ±5.6 


15% 


30 


78.05 


1.29 


(1) 


-49.12 ±0.14 


-160.46 ±0.11 


-6.5 ±9.1 


-14.4 ± 10.9 


67% 


31 


72.37 


1.09 


(1) 


-62.09 ± 0.09 


-76.67 ±0.11 








32 


80.37 


0.76 


(4) 


-57.43 ± 0.09 


-162.26 ±0.11 








33 


79.53 


0.75 


(3) 


-32.44 ±0.10 


-29.99 ±0.12 


-4.5 ±3.4 


-12.4 ± 3.7 


28% 


34 


80.79 


0.71 


(1) 


425.21 ±0.10 


158.23 ±0.11 








35 


72.16 


0.67 


(1) 


-60.58 ±0.10 


-77.28 ±0.11 








36 


78.48 


0.59 


(2) 


-49.46 ± 0.08 


-161.70±0.13 








37 


66.89 


0.56 


(1) 


-66.88 ±0.12 


-61.05 ±0.11 








38 


80.58 


0.55 


(1) 


437.41 ±0.10 


67.80 ±0.15 


11.4 ±7.2 


-6.8 ± 10.8 


63% 


39 


75.32 


0.49 


(2) 


-548.55 ± 0.09 


-314.11 ±0.11 








40 


79.53 


0.48 


(2) 


-49.38 ± 0.09 


-162.86 ±0.17 








41 


76.79 


0.44 


(2) 


-548.30 ± 0.09 


-314.13 ±0.09 








42 


79.32 


0.43 


(4) 


-293.32 ±0.08 


-100.02 ±0.08 








43 


67.94 


0.42 


(1) 


-53.79 ±0.09 


-54.67 ±0.11 








44 


80.58 


0.39 


(1) 


-287.44 ±0.10 


-98.35 ±0.12 


-18.6 ±3.5 


-1.5 ±3.8 


19% 


45 


88.59 


0.38 


(1) 


-64.55 ±0.13 


-162.29 ±0.12 


-7.6 ±8.2 


-34.9 ± 10.7 


30% 


46 


63.94 


0.38 


(4) 


-54.44 ±0.11 


-54.30 ± 0.09 








47 


80.16 


0.36 


(1) 


440.10 ±0.10 


58.11 ±0.16 


18.7 ±7.1 


-26.3 ± 10.6 


30% 


48 


62.47 


0.35 


(4) 


-64.15 ±0.09 


-60.07 ± 0.09 








49 


79.32 


0.33 


(3) 


-293.50 ±0.10 


-100.37 ±0.12 


-8.8 ±3.5 


-8.5 ± 4.0 


31% 


50 


78.26 


0.32 


(2) 


-61.72 ±0.08 


-78.38 ±0.17 








51 


62.25 


0.31 


(3) 


-54.39 ±0.10 


-54.57 ±0.10 








52 


69.42 


0.31 


(1) 


-40.40 ±0.10 


-7.34 ±0.11 








53 


79.74 


0.31 


(1) 


284.39 ±0.10 


-55.26 ±0.12 


0.2 ±4.0 


-18.3 ±4.9 


27% 
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Table 2. continued. 



Feature 


Vlsr 
(kms 


Fpeak 

(Jy beam"') 


Ax 
(mas) 


Ay 
(mas) 


(kms ') 


(kms 


A 1 T 71 / 1 T 71 

A|y|/|V| 


54 


77.21 


0.30 (1) 


-58.94 + 0.10 


-77.13 ±0.13 








55 


80.16 


0.28 (1) 


437.51 + 0.10 


67.22 ± 0.20 








56 


76.16 


0.27 (1) 


-103.16 ±0.10 


-150.22 ±0.12 








57 


89.43 


0.26 (1) 


-49.16 + 0.12 


-158.89 ±0.12 








58 


71.95 


0.23 (1) 


-67.77 ±0.10 


-61.37 ±0.11 








59 


80.58 


0.22 (4) 


-82.16 ±0.08 


-104.88 ±0.09 








60 


77.42 


0.21 (2) 


-546.43 ± 0.08 


-314.43 ±0.10 








61 


70.68 


0.21 (2) 


-65.69 ± 0.08 


-60.78 ± 0.09 








62 


77.21 


0.20 (1) 


418.75 ±0.10 


92.99 ±0.12 


5.8 ±3.8 


-1.6 ±4.7 


64% 


63 


70.05 


0.20 (3) 


27.36 ± 0.09 


-15.57 ±0.12 








64 


59.51 


0.19 (3) 


-57.57 ± 0.09 


-63.38 ±0.11 








65 


75.74 


0.17 (3) 


-549.28 ± 0.09 


-3 14.97 ±0.11 








66 


81.63 


0.16 (1) 


-70.96 ±0.10 


-88.63 ±0.13 








67 


79.74 


0.16 (3) 


-53.11 ±0.09 


-105.70 ±0.10 








68 


70.89 


0.16 (4) 


-51.85 ±0.09 


-73.26 ±0.10 








69 


74.68 


0.15 (4) 


31.19 ±0.09 


-13.67 ±0.09 








70 


70.89 


0.14 (4) 


-38.58 ± 0.09 


-53.28 ± 0.09 








71 


68.57 


0.14 (4) 


-53.90 ±0.09 


-54.56 ± 0.09 








72 


67.31 


0.13 (1) 


-63.71 ±0.10 


-60.09 ±0.14 








73 


75.10 


0.13 (1) 


442.20 ±0.10 


98.34 ±0.14 








74 


75.74 


0.12 (1) 


-51.19±0.10 


-68.44 ±0.16 








75 


66.05 


0.11 (3) 


-541.56 ±0.10 


-291.69 ±0.12 








76 


79.53 


0.11 (4) 


-51.08 ±0.09 


-160.20 ±0.10 








77 


73.21 


0.10 (1) 


-67.84 ±0.10 


-61.58 ±0.15 








78 


81.00 


0.09 (4) 


412.27 ± 0.09 


171.84 ±0.10 








79 


73.84 


0.09 (2) 


-68.81 ± 0.09 


-61.71 ±0.15 








80 


76.16 


0.09 (3) 


c A r\ A r\ r\c\ 

-540.64 ± 0.09 


-261 .93 ± 0.14 








81 


77.21 


0.09 (1) 


419.08 ±0.10 


168.99 ±0.13 








82 


71.94 


0.09 (3) 


-99.30 ± 0.09 


-173.21 ±0.11 








83 


69.42 


0.09 (1) 


-68.13 ±0.11 


-61.87 ±0.15 


-26.8 ± 8.2 


8.0 ± 10.6 


30% 


84 


76.58 


0.08 (1) 


-52.05 ±0.10 


-68.04 ±0.16 








85 


63.73 


0.08 (2) 


-53.41 ±0.09 


-54.35 ±0.14 








86 


85.85 


0.07 (4) 


-4.96 ± 0.09 


-151.45 ±0.12 









Reference Feature: Absolute position & Proper motion 



Feature 


R.A.(J2000) 


Dec.(J2000) 


Vx 




A|V|/|V| 




(hm s) 


(°' ") 


(kms-i) 


(kms-') 




24 


18:34:40.28950±0.00010 


-09:00:38.2044±0.0015 


5.3 ±8.6 


21.3 ±8.6 


39% 
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Table 3. Parameters ofVLBA 1.665 GHz hydroxyl maser features. For each identified feature, the label (given in Col. 1) increases 
with decreasing brightness. Cols. 2(3) and 4(5) report the LSR velocity and brightness of the brightest spot for the right (left) circular 
polarization. Cols. 6 and 7 give the position offset relative to the feature # 1, toward the East and the North directions, respectively. 
Cols. 8 and 9 report the fit split of LSR velocity (RCP-LCP) between co-spatial circular polarization components and the inferred 
1.0. s. magnetic field strength, respectively. The absolute position of the reference feature # 1 is reported at the bottom of the Table. 



Feature 


(kms^') 
RCP LCP 


Fpeak 

(Jy beam~') 
RCP LCP 


Ax 
(mas) 


Ay 
(mas) 


Ay/ 

(km s"') 


Bz 
(mG) 


1 


70.88 


71.76 


2.39 


0.40 








-0.86 


-1.5 


2 


67.89 




2.02 




-210.26 ± 1.07 


178.74 + 0.77 






3 


66.84 


66.84 


1.36 


0.84 


14.24 ± 0.62 


18 87 + 77 






4 


73.87 


73.87 


1.01 


0.16 


-251 52 + 78 


88 39 + 1 47 






5 




74.75 




0.97 


-385 08 + 1 46 


203.14 ± 1.00 






6 




70.57 




0.89 


199.12 ± 1.45 


174 08 + 1 13 






7 


68.07 




0.74 




-185 97 + 67 


191 75 + 1 29 






g 


72.99 


72.99 


0.67 


0.26 


57.29 ± 0.71 


119.21 + 0.94 






9 




75.10 




0.60 


-297.88 ± 1.46 


130.40 + 1.14 






10 


74.92 




0.55 




-358 64 + 1 31 


-170.67 + 1.14 






11 


70.35 


69.82 


0.52 


0.23 


56.12 ± 0.94 


114.09 + 1.10 


+0.62 


+1.0 


12 


67.19 




0.46 




4-; QO + 66 


111 93 + 08 






13 


76.33 


79.67 


0.25 


0.43 


-223.40 + 1.45 


71.08 + 1.20 


-3.41 


-5.8 


14 




76.33 




0.41 


366.46 ± 1.57 


-149.68 ± 1.13 






15 


68.77 




0.38 




52.32 ±2.18 


113.99 + 1.27 






16 




63.50 




0.31 


32.32 ± 1.56 


108.87 + 1.73 






17 


64.90 




0.30 




39.23 ± 1.10 


111.08 + 1.12 






18 


70.18 




0.25 




192.37 ± 1.33 


175.05 + 1.53 






19 


79.32 


77.21 


0.25 


0.17 


-321.93 + 0.90 


130.00 + 2.23 


+2.11 


+3.6 


20 


64.55 


65.08 


0.23 


0.13 


37.16 + 0.92 


108.63 + 1.33 


-0.49 


-0.8 


21 


66.84 




0.22 




-134.66 + 0.91 


147.73 + 1.29 






22 


62.62 


62.79 


0.16 


0.20 


26.50+ 1.98 


106.74 + 1.81 






23 


69.65 


69.30 


0.18 


0.19 


17.17 + 1.78 


18.23 + 2.17 






24 


60.68 


60.51 


0.12 


0.14 


20.64 ± 1.70 


70.44 + 2.47 







Reference Feature: Absolute position 



Feature 


R.A.(J2000) 


Dec.(J2000) 




(hm s) 


("' ") 


1 


18:34:40.28515+0.00026 


-09:00:38.3783+0.0044 



(a) A magnetic field of 0.34 mG shifts the emission velocity of 1665 MHz OH masers by 0.2 km s ' 
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Table 4. Parameters of EVN 6.7 GHz methanol maser features. For each identified feature, the label (given in Col. I) increases with decreasing 
brightness. Cols. 2, 3 and 4 report the LSR velocity and brightness of the brightest spot, and its percent brightness variability. Cols. 5 and 6 give the 
position offset relative to the feature # 10, toward the East and the North directions, respectively. Cols. 7 and 8 report the projected components of 
the feature proper motion relative to the center of motion (as defined in Sect. l4.4l and identified with label # 0), along the East and North directions, 
respectively. Col. 9 shows the relative error of the proper motion amplitude. "S" is the barycenter of the methanol maser distribution as defined in 
Sect. 14.41 The absolute position of the reference feature # 10 is reported at the bottom of the Table. 



Feature 


TT 

'LSR 


P 


Var." 


Ax 


Av 


' X 




AIVI/IVI 




(km s"^) 






( mil 


( TTlfl 


(km s"') 


(km s"') 




1 


74 79 

/t. / Z 


1 4f^ 0^ 




1 44 SO + 08 

iTH-.JW IE \J.\JO 


^ 1 q 1 + n OQ 


1 Q + Q 
i ic \j.y 


/I 8 + Q 
T-.o ic yj.y 


1 Q% 

ly /c 


9 


7S SI 

/ J.J 1 


1 1 ^ 88 


71 f?^ 


1 77 1 7 -1- 07 
i / / . 1 / It w.w / 


7^^ Ql + 08 








J 


74 98 


10^ 74 


1 \J W 


"^Afx 87 + 08 


zl7S QO + n 08 

T-ZJ.7W 3: W.WO 


^ A + Q 
j.o ic \j.y 


— 1 7 + Q 

— i .z ic yj.y 


94% 

Zt- /O 


A. 


74 QS 


1 07 78 


LA/0 


1 8^S 40 + 08 

IoO.'tW It VJ.Uo 


7^7 + 08 
ZJ / .y J x W.WO 


7+08 
W. / IC w.o 


7 7 + Q 

/ . / IC W.7 


1 7% 

i z /o 




74 81 

/ T.O 1 


80 f\A 


Jo /O 


1 f^8 1 7 -1- 08 

iOO. i / It U.Uo 


7/Ll 8Q + OQ 








D 


80 f^Q 


J J.ZO 


fi7% 


—47 Q8 -1- 08 
— '-r 1 .yo It U.Uo 


zL3 87 + 08 

H-J.O / 3: W.WO 


—74 + 08 
— z.T- IC w.o 


8 7 + Q 
o.z ic yj.y 


1 1 % 

i i /o 


7 


7S 8f^ 

/ J.OU 


47 A'~K 


7zL% 


1 7f^ 77 + 08 
i / U.Z / It u.uo 


7^Q OS + 08 

ZJ7.WJ 3: W.WO 








Q 
O 


81 74 


^7 47 


1 8% 
i /o 


1 70 -1- 07 
i .ZW It w.u / 


1 77 + 08 

i .Z / IC W.WO 


—34 + 08 
— J .T- IC W.O 


—0 Q + 8 
— yj.y ic w.o 


73% 
ZJ /c 


Q 


74 07 

/ T. WZ 






1 Sf^ 8S + 08 
i JU.oJ It W.Wo 


30/L 3S + OQ 
JWt-.JJ x W.W7 


zl 7 + Q 
H-.z ic \j.y 


4 3 + 10 
J IC i .w 


1 fi% 

io /o 


in 


87 "Kfs 


IQ 








w 


-S 7 + 8 
J . / IC w.o 


—07+08 
w.z IC w.o 


14% 

iT- /C 


1 1 


79 70 


7f^ 


71 


77^ 1 S + 08 
Z / J . 1 J It U.Uo 


777 3zL + 08 

Z / / . JH- IC W.WO 


T 7 + 8 
J.Z ic w.o 


— 3 8 + Q 
— J.o ic yj.y 


1 8% 
i /o 




81 91 
1 .z 1 


^4 7^^ 


zo /o 


_'17 71 + 07 

— J / .Zl ItW.W/ 


f,^ 77 + 08 

i . / / x W.WO 


— T A + 8 
— J .0 ic w.o 


— T 1 + Q 
— j.i ic yj.y 


1 8% 
i /o 




74 4S 


9S 8f^ 




9 1 ^ 7Q _i_ n OQ 

Z i J . / !7 It W.W" 


778 9S + 1 1 

Z/O.ZJ IC W.il 








It- 


74 78 


99 88 
zz. 00 




970 S8 + 08 
ZZW.JO It w.wo 


77s 38 + 10 
Z/J.JO 3: W.iW 








1 ^ 


74 


91 S4 


J\J/C 


^OQ 47 + 17 

JW^^.H-Z It U.iZ 


zl7S 1 fi + 1 fi 


1 S Q + 1 f, 
ij.y ic i .0 


_1 A + 1 Q 
i .0 ic L .y 


1 0% 
iw /o 


\fx 
lO 


7Q 4^1 


70 f\\ 




—^4 Q8 + 08 

— Jt.!70 It U.Uo 


Al Ad + D 08 

H / .T-H IC W.WO 


— T 7 + Q 
— J .z ic yj.y 


S 4 + Q 
j.T- ic yj.y 


1 S% 

ij /o 


17 
1 / 


1^ 84 


1 S 77 


1 10 


181 Q8 + 08 
101 ."0 It U.Wo 


777 7zL + OQ 

1 Z^. 1 ^ IC W.W" 


7 + Q 
w. / ic yj.y 


— 1 S + 1 

— i . J IC i .w 


ow /c 


1 R 


7^ 7S 


1 ^ S8 




744 f\\ + OQ 
ztt.ui It \j .yjy 


77A 71 +0 10 

Z./\J.Z.L ICW.iW 


17 1+10 

i Z . i IC i .w 


3 7 + 17 

J .Z IC i .z 


Q% 
y /c 


1 Q 


7^ 4Q 


1 ^ 4^ 


Jo /o 


\fn 1 4 + 08 

iU / . It It U.UO 


781 S7 + OQ 
zo i . J / ic yj.yjy 


— 7 A + Q 
— z.o ic \j.y 


—7 1+10 
— / . i IC i .w 


1 3% 

i J /o 


90 

Z-U 


74 IQ 


1 7 QO 


uw /c 


1 QzL 7 1 + OQ 
ly-t.Al It \J.\Jy 


111 zLQ + 13 

Z/Z.T-7 IC W.iJ 








91 


7f=i 74 

/ 1). /t 


19 77 
iz. / / 




1 QO 9^ -1- 07 
17U.ZJ It yj.yj/ 


77Q QQ + 08 
ZAy .yy ic w.wo 


10 + 08 
iW.W ic w.o 


_1 7 + Q 
i .z ic yj.y 


/o 


99 


74 8Q 


19 98 
1 z.zo 


73% 

f J /c 


'\AA QS + 1 Q 
j'-r'-r.yj It yj.ly 


zLl 8 87 + 71 

T-iO. OZ IC W.Zi 


14 fi + 7 7 

iH-.U IC Z. / 


— S Q + 7 7 
— J .y rji z. 1 


1 7% 

i / /o 


9^ 

Zj 


74 ^7 


1 1 81 


00 /o 


70S QO + 17 

ZWJ."U It W.iZ 


7fi7 3fi + 1 Q 
zo / .JO IC W. i 7 








94 


74 09 

/ T. wz 


1 (\f\ 


" /c 


9^A S9 + OQ 
ZJU.JZ It \J.\Jy 


78Q 8fi + 1 1 

ZO7.OO IC W. i i 








9^ 


80 f^O 


1 01 




^^1 SQ + 08 

JJi.J!7 It W.WO 


— 1 zl3 zl3 + OQ 

— IH-J.H-J IC W.W7 


4 3 + 10 
^. J IC i .w 


—f, 8 + 11 
— 0.0 IC i . i 


1 4% 

iH- /O 




7^ 9^ 

/ J.Z J 


Q 7S 




7^ SS + 08 
/ J .J J It w.wo 


3zL1 OS + OQ 
Jt-1 .WJ IC w.w^ 


— ^ 1 + Q 
— J . 1 ic yj.y 


7 0+10 
/ .W IC i .w 


1 3% 
i J /o 


97 


7Q 98 


8 7^ 


3Q% 


7AS 6^ + 08 

ZtJ.UJ 31 W.WO 


-S4. 11 + OQ 

JT^. i i IC W.W" 


T + Q 
J .w IC \j.y 


— S S + 1 

J — J IC i .w 


i U /CJ 


9R 
zo 


81 S7 


8 IQ 


1 3% 


— 4Q 08 + 08 
— '-ty .\Jo It U.Uo 


1 03 30 + OQ 
iWJ.JW IC w.w^ 


— A Q + Q 
— \j.y ic yj.y 


7 4+10 

Z.T- IC i .w 


1 3% 

i J /o 


9Q 

Z7 


7^^ 9^ 


f\ 81 

U.O i 


78% 


770 f^7 + 08 

ZZW.UZ It W.WO 


77^ Q7 + 10 

ZZJ.^J IC W.iW 


3 3 + 10 
J.J IC i .w 


— 17 + 1 1 

— i .Z IC i . i 


7Q% 

Z7 /c 


D\J 


74 79 

/ T. / Z, 


6 S7 

U.J / 




1 Q 1 Q4 + 7 1 
ly 1 .y-r X w.z 1 


77Q 1^ + 07? 

Z/I/.iJ IC w.zz 








'^1 


7^ 9^ 


9^ 

U.Z J 




^1 S 81 -1- 08 

J 1 J .0 1 X W.WO 


37S zLS + OQ 
JZJ.H-J IC W.W7 


13 + 10 
i . J IC i .w 


—70+1 9 

— Z.W IC i .z 


48% 

T-O /o 


JZ 


lA 1 

/ T. i W 


S 8f^ 

J.OU 


3S% 

JJ /c 


99Q 71 -1-0 11 

ZZO . / i It W. i i 


78/1 8/1 + 1 

ZOt-.OH IC W.iO 










80 8f^ 


S 1 S 
J. i J 


1 01 % 


A SQ + 1 1 

'-r.Oy It W. i i 


8 87 + 1 S 
0.0 / IC W. i J 








Ot 


77 7^^ 
/ / .zu 


4 

T-. J J 


31 % 


1 Qf^ Sf^ + 08 

i!7U.JU It W.WO 


731 47 + OQ 

Z J i .H-Z IC W.W" 










74 78 


4 

T.J J 




1 Q 1 OQ + 14 

1 " 1 .W" It W. It 


78Q 30 + 7S 

ZO7.JW IC W.ZJ 








D\J 


74 10 

/ T. i \J 


4 ^0 

T. JW 




1 + 11 

1 J ."J X W. 1 1 


7S8 Q3 + Ifi 
zjo.yj IC W.iU 








/ 


87 f^9 
oz.uz 


^ 7^ 

J.Z J 


1 S% 


78 -1- 08 

— JJ. / It W.WO 


/lO 9Q + OQ 
T-W.Z7 IC W.W7 


—4 3 + 10 
— ^. J IC i .w 


S 9 + 1 7 

J .Z IC i .z 


1 7% 
i / /c 


'^S 
Oo 


8^ 9^ 
J.Z J 


^ 1 8 
J. i 


38% 
JO /o 


— 4Q 41 + 08 

— T!7.Ti It W.WO 


Q7 Q^ + OQ 

" / .7 J IC W.W7 


— S S + Q 
— J .J ic yj.y 


13 + 10 
i . J IC i .w 


1 7% 

i / /o 




80 9S 


7 81 


3Q% 


^1 ^ 77 + 08 

JIJ.Z,/ 31 W.WO 


_1 7S SQ + OQ 

i JJ.J" IC W.W" 


1+10 

W. i IC i .w 


— A 7 + 1 7 

U. / IC i .z 


1 8% 
i /c 


tU 


7^ ^1 


7 f^S 

Z.O J 


1 3S% 

i J J /C 


1 8S ^8 -1- OQ 

lOJ.JO It W.W" 


70Q 7S + 1 1 
z,\jy . / J ic w. i 1 








41 

Ti 


7^ 


7 S8 

Z.Jo 




1 7^ ^0 -1- 1 ^ 

1 / U. JW It W. 1 J 


7fiO 1 3 + 7Q 
ZOW.iJ IC W.Z7 








42 


72.79 


2.27 


29% 


262.25 + 0.10 


221.26 + 0.15 








43 


81.48 


2.06 


77% 


-51.58 + 0.09 


39.14 + 0.12 








44 


78.93 


1.94 


24% 


78.24 + 0.08 


-52.01+0.09 


-4.2+1.0 


1.4+1.2 


24% 


45 


73.58 


1.85 




252.27 + 0.11 


281.21+0.19 








46 


81.04 


1.79 


11% 


-45.07 + 0.09 


-31.48 + 0.12 


-8.1 ± 1.1 


0.1 + 1.5 


15% 


47 


78.84 


1.78 


124% 


-145.32 + 0.15 


294.05 + 0.23 


-5.7 + 1.8 


11.4 + 2.5 


19% 


48 


73.23 


1.72 


136% 


117.85 + 0.26 


315.86 + 0.44 








49 


72.87 


1.71 


28% 


331.36 + 0.12 


326.56 + 0.17 


2.0 ± 2.0 


-5.4 + 2.5 


42% 


50 


80.42 


1.65 


67% 


118.38 + 0.09 


-15.13 + 0.12 


-3.5 + 1.7 


-2.0 + 2.1 


45% 


51 


72.17 


1.29 


21% 


189.08 + 0.09 


215.22 + 0.11 


1.2 + 1.2 


-1.6 + 1.4 


65% 


52 


80.95 


1.18 




285.89 + 0.16 


-92.10 + 0.17 








53 


78.58 


1.08 


76% 


361.62 + 0.16 


-60.32 + 0.21 
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Table 4. continued. 



Feature Vlsr 


Fpeak Var.« 


Ax 


Ay 


Vx 


Vy A\V\/\V\ 


(kms-i) 


(Jy beam"') 


(mas) 


(mas) 


(km s-i) 


(km s-i) 



54 


79.90 


1.04 


58% 


55 


70.24 


0.91 


14% 


56 


83.06 


0.88 


25% 


57 


73.58 


0.84 




58 


82.44 


0.83 


18% 


59 


78.49 


0.78 


32% 


60 


69.80 


0.76 


45% 


61 


76.56 


0.69 




62 


79.81 


0.60 


57% 


63 


80.16 


0.55 


40% 


64 


80.16 


0.54 


28% 


65 


78.76 


0.47 


7% 


66 


76.91 


0.47 




67 


79.72 


0.46 


42% 


68 


78.23 


0.41 


18% 


69 


80.25 


0.37 


62% 


70 


77.44 


0.36 


65% 


71 


77.70 


0.29 


85% 


72 


78.76 


0.25 


20% 


73 


78.23 


0.25 


23% 


74 


79.20 


0.24 


46% 


75 


77.88 


0.22 


45% 


76 


78.14 


0.19 


40% 


77 


78.14 


0.17 


54% 


78 


77.35 


0.16 


47% 


79 


72.35 


0.13 


28% 


80 


78.32 


0.12 


5% 


81 


67.87 


0.10 









148 


46 


+ 





11 


-102 


60 


+ 





16 


192 


80 


+ 





09 


47 


85 


+ 





11 


202 


89 


± 





10 


-1 


49 


± 





13 


221 


83 


+ 





13 


263 


46 


+ 





18 


211 


57 


+ 





25 


-10 


52 


+ 





33 


365 


62 


+ 





10 


-53 


08 


+ 





16 


-54 


83 


± 





11 


303 


27 


± 





15 


84 


27 


± 





24 


381 


13 


± 





26 


110 


44 


+ 





16 


-67 


92 


+ 





24 


109 


85 


+ 





18 


-72 


17 


+ 





26 


369 


64 


+ 





19 


-142 


39 


+ 





24 


18 


53 


± 





14 


298 


80 


± 





22 


165 


18 


± 





25 


-23 


41 


± 





28 


160 


89 


+ 





16 


-26 


57 


+ 





27 


83 


31 


+ 





13 


-38 


89 


+ 





18 


354 


72 


+ 





21 


-140 


73 


+ 





38 


-30 


76 


± 





14 


69 


63 


± 





22 


-142 


36 


+ 





26 


302 


09 


± 





41 


1 

— i 




+ 


yj 




1 

jLj 


'20 
JO 


+ 


u 


'20 
JO 


167 


15 


± 





19 


-25 


88 


+ 





33 


-65 


58 


+ 





19 


102 


70 


+ 





32 


-69 


95 


± 





33 


103 


10 


± 





44 


200 


35 


+ 





35 


178 


88 


+ 





46 


-68 


92 


+ 





37 


103 


28 


+ 





62 


-151 


21 


± 





28 


303 


05 


+ 





43 


178 


61 


± 





38 


159 


23 


± 





56 


50 


77 


+ 





32 


257 


35 


+ 





50 


178 


94 


+ 





48 


49 


84 


+ 





55 


131 


30 


± 





07 


149 


69 


± 





08 



155 143 



-3.3 + 1.5 
-2.4 ± 1.3 
1.5 ± 1.7 


-11.1 ±2.0 
-5.2 ± 1.6 
-3.7 ±2.1 


17% 
28% 
51% 


-5.4 ± 1.5 


10.0 ± 2.0 


17% 


-2.2 ±3.0 
1.4 ±2.4 


-10.4 ±3.2 
7.4 ± 3.2 


30% 
43% 


-1.9 ±2.3 


-0.8 ± 2.8 


114% 


-9.4 ± 3.7 


9.1 ±5.1 


34% 







Reference Feature: Absolute position 



Feature 


R.A.(J2000) 


Dec.(J2000) 




(hms) 


r ") 


10 


18:34:40.27658±0.00014 


-09:00:38.4086±0.0021 



(a) VariabiUty (%) = (Fmax - Fmin of the brightest spot) / ((Fmax + Fmin)/2.0) 



